Electronic Structure of Calcium Hexaboride within the Weighted 

Density Approximation 

Zhigang Wu 1 , D.J. Singh 2 and R.E. Cohen 1 
1 Geophysical Laboratory, Carnegie Institution of Washington, 
5251 Broad Branch Road, NW, Washington, DC 20015 
2 Center for Computational Materials Science, 
Naval Research Laboratory, Washington, DC 20375 
(Dated: February 2, 2008) 

Abstract 

We report calculations of the electronic structure of CaBg using the weighted density approx- 
imation (WDA) to density functional theory. We find a semiconducting band structure with a 
sizable gap, in contrast to local density approximation (LDA) results, but in accord with recent 
experimental data. In particular, we find an X-point band gap of 0.8 eV. The WDA correction of 
the LDA error in describing the electronic structure of CaBg is discussed in terms of the orbital 
character of the bands and the better cancelation of self-interactions within the WDA. 
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The discovery of ferromagnetism in the lightly doped CaBg and SrBg resulted in con- 
siderable experimental and theoretical work aimed at understanding this unexpected phe- 
nomenon. The reported ferromagnetism is particularly remarkable considering that the 
material, Ca 1 _ :r La a: B 6 , contains no elements with partially filled d or / orbitals, and even 
more strangely, the ferromagnetism is associated with very low doping levels (< 1%), and 
therefore low carrier concentrations, has a very low moment (~ 0.07 /ie/La for x=0.005), 
but a very high Curie temperature, Tc as high as 600K (Ref. |l|). Further, a semiconducting 
material, with ferromagnetism of this type, could potentially be of importance in producing 
room temperature spintronic devices. [2J, y| 

Two different classes of theories were advanced early on in an effort to explain the strange 
ferromagnetism. The first was based on the doping of an excitonic or like system, i.e. an 
unconventional state in the parent, CaBg, J, 0j6,_Z| and the second was based on the 
ferromagnetic phase of the dilute electron gas. l|, lalSl In the first class, the ferromagnetism 
arises from La induced doping of carriers into an excitonic ground state. This depends on a 
bare band structure that is semimetallic or very close to it, so that the prerequisite instability 
of CaBg to an excitonic state can exist. The second class mentioned does not depend 
on such an electronic structure, but rather involves a ferromagnetic instability of a low 
density electron gas, though this class of theories suffers from difficulties obtaining the very 
high reported values of Tq. More pedestrian explanations for the observed ferromagnetism 
involving defects and their interactions and also extrinsic effects have been advanced more 
recently. 



This naturally resulted in a strong interest in the electronic structure of the parent com- 
pound CaB 6 . First-principles calculations within the local density approximation (LDA) 
predicted that CaBg and SrBg are semimetals with a small overlap at X point of the sim- 
ple cubic Brillouin zone, or very small gap insulators, jl^l I5I ] consistent with the excitonic 
theories being formulated at that time. However, it is well known that the LDA often 
underestimates band gaps, and in some extreme cases it can predict a semimetallic band 
structure for a semiconductor, one example being Ge. [16 ! ] In simple semiconductors, the 
LDA underestimates of the band gaps are understood as being associated with an extreme 
non-analytic and non-local behavior of density functional theory as the particle number is 
changed to include a single conduction electron, \v\ . implying the need for a self-energy 
correction, or at least an orbital dependent potential to obtain a realistic description of band 
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gaps. 

The many-body quasiparticle GW approximation [19( is very succesful in describing the 
band structures of insulators. In particular the self-energy correction implicit in it corrects 
most of the band gap errors in simple semiconductors and other materials as well. 2o| 
However, its results for CaB 6 have differed depending on the formulation of GW being used. 
Pseudopotential based GW calculations j| yielded a sizable 0.8 eV band gap at X point, 
while recent all-electron GW results showed an unexpected increase in band overlap 



relative to the LDA, perhaps related to the treatment of deep-core states. |22j, Ul 

The electronic structure of CaBg has also been the subject of some debate from an ex- 
perimental point of view, although a semiconducting electronic structure has become more 



enerally accepted now. de Haas-van Alphen (dHvA) [2J, 25] and Shubnikov de Hass (SdH) 
2y| experiments were consistent with a semimetallic nature, while recent angle-resolved 
photoemission spectroscopy (ARPES) and resonant inelastic x-ray scattering (RIXS]mea- 
surements imply that CaP>6 is a semiconductor with a 1.0 (Ref. l27l ) - 1.15 (Ref. 1281 1 eV 
band gap. Recently, Rhyee et al. reported transport, optical, and tunneling measurements 
all consistent with a 1.0 eV band gap for pure CaBg. They argued that the previously mea- 
sured semimetallic characteristics in CaB 6 , SrB 6 , and EuB 6 actually originate from boron 
vacancies. 29f 

We note that, as pointed out by Rodriguez and co-workers, the valence and 

conduction bands of CaBg and the related hexaborides are of very different character, namely 
B p valence bands and Ca derived conduction bands. This is reminiscent of the situation 
in the rare-earth trihydrides, e.g. YH3. Those materials feature a semiconducting ground 
state, which is however predicted to be metallic in the LDA due to a band overlap between 
nominally occupied H derived states and nominally unoccupied metal states. In YH3, this 
error is understood to be due to incomplete cancelation of self-interaction errors on the H 
sites in the LDA. These self-interactions lead to a H band position that is too high relative 



to the conduction bands, which is at least partly corrected in GW calculations. 
It should be noted though that such self- interactions are a local effect, more akin in some 
ways to the on-site Coulomb effects that lead to Mott insulating states in transition metal 
oxides, than to the non-local, non-analytic discontinuities in density functional theory that 
are important in simple semiconductors. 3^ As such, the rare earth trihydrides and perhaps 
CaBg are interesting test cases for theories of correlated electron systems, including density 



3 



functional based approaches. 

Here, we present results of first-principles calculations based on the density functional 
theory (DFT) within a non-local density approximation, namely the weighted density ap- 
proximation (WDA). 35, 36, 37, |^ The WDA uses a model pair-distribution function to 
incorporate non-local information about the charge density into the exchange correlation 
energy and potential, via an exact expression for the total energy. Unlike GW and like 
approaches, the WDA is oriented towards total energies rather than excitation spectra. As 
such, not surprisingly, it improves ground state properties, like lattice parameters, relative 
to the LDA. Like the LDA, the WDA is exact for the total energy of the uniform electron 
gas, but unlike the LDA it is also exact for single electron systems, where it is self-interaction 
free. We emphasize that the WDA is a smooth non-pathalogical approximation to density 
functional theory, with only short range non-locality. The WDA exchange correlation poten- 
tial is a standard orbital independent (for valence states) potential, which different orbitals 
feel differently only through differences in the spatial distrubutions of their charge densities. 
It does not include the discontinuities that are present in the exact density functional theory 
and are responsible for most of the band gap correction in simple semiconductors like Si. 
Detailed calculations for YH 3 and LaH 3 have shown that the WDA can yield insulating 
gaps of similar magnitude to the GW approximation in these hydrides due to the better 
cancelation of self-interaction errors, while at the same time improving the quality of ground 
state properties. j^J 

In the WDA the exchange-correlation (xc) energy expression is very similar to the general 
DFT form, 



ET>} =j j ^^G[|r-r'U(r)]^r', (1) 

where G is the approximate model coupling constant averaged pair-distribution function, 
and n(r) is the weighted density which can be determined from the sum rule which states 
that the charge density in the exchange correlation hole is unity. 

J n{r')G[\r-v'\,n{r)\dr' =-l. (2) 

The WDA is exact for both a one-electron system and the uniform electron gas, independent 
of the model G used, provided that it is monotonic and properly normalized. Because the 
exact G function is not known, different forms of G have been tried, 40] from which we 
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chose the following one 

G(r,n) = c{l-exp(-[^]- 4 )}, (3) 
This choice yields good energetics (structural properties) and band structures for the rare 



earth trihydrides. 



39| This type of G was derived from the Gunnarsson- Jones 4l( ansatz, 



but it has a shorter range than the orig inal one. Details of the WDA formalism used in this 



paper can be found in Refs. 



39 and 



42 



The calculations were done using a planewave basis with pseudopotentials. The k-grids for 
the Brillouin zone samplings during iteration to self-consistency were as dense as 12 x 12 x 12, 
and the energy cutoff was 121 Ry. 43] The Ca semicore states were treated as valence states 
for the purpose of constructin g T roullier-Martins pseudopotentials. Shell partitioning was 
employed as discussed in Ref. 38. The crystallographic structure of CaB 6 is simple cubic, 
with Ca atoms located at cubic corners. The B atoms form octahedral cages around the 
cube center. The internal parameter x determines the distance between a B atom and the 
cubic center. When x = 0.207, the intra- and inter- octahedral B-B distances are equal. 
Our calculated lattice constant a and internal parameter x within the LDA and WDA are 
shown in Table |U The WDA results are in excellent agreement with experimental data, 
while the LDA predicts a lattice constant 1.6% smaller, which is a typical LDA tendency. 
Thus, reflecting its origins in density functional theory, where energy and density are primary 
variables, the WDA improves upon the LDA structure. 

To facilitate comparison of the LDA and WDA electronic band structures, we used the 
experimental values of a and x for both the LDA and WDA. Fig. ^ shows our results 
for the band structures of CaBg. The LDA (left panel) predicts a tiny band gap (slightly 
more than 0.1 eV) at t 



he X-point, which agrees with Ref. j2JJ, while other LDA results 
show a tiny overlap llfil l. This small difference in band structure is due to different 

choices of the LDA scheme, chemistry (Sr vs. Ca) and lattice structure as well as the 
approximations in using pseudopotentials. To check this we did well converged general 
potential linearized augmented planewave (LAPW) calculations within the LDA for two 
structures. The first is the same experimental structure as used for Fig. and the second 
is a structure where we relaxed the internal coordinate x in the LDA keeping the lattice 
parameter fixed at the experimental value. This yielded x=0.2025. With the experimental 
value of x=0.207 the LAPW calculations yielded a small band gap of 0.2 eV, while at 
x=0.2025 a semimetal with a 0.1 eV band overlap is predicted. This implies that the 



structure is the key ingredient in reconciling the various results and also gives an error 
estimate of ~ 0.1 eV for our pseudopotential calculations. On the other hand, the right 
panel shows that within the WDA, the bands are shifted, quite similar to the cases of YH 3 
and LaH 3 , resulting in a sizable band gap of 0.8 eV at X point. Parallel WDA calculations 
using the uniform electron gas G function also yield a gap opening, though by a smaller 
amount. The band gap in that case is 0.3 eV. 

Thus the WDA, with an appropriate choice of G, predicts CaB 6 to be a semiconductor 
with a gap comparable to that found in standard GW calculations [3] and in experimental 
spectroscopic measurements. 22, 2slQ| This shows that there exists a physical approxima- 
tion to density functional theory that is continuous and has only short range non-locality 
and that gives both a reasonable band gap and ground state properties (crystal structure) in 
agreement with experiment for CaBg. We expect that WDA calculations would yield band 
gap opening in SrB 6 as well. 

In summary, WDA calculations predict that CaBg is a semiconductor, and the calculated 
band gap is in excellent agreement with recent experiments and quasi-particle GW calcula- 
tions. Thus the WDA improves both the structural properties and electronic structure of 
these hexaborides, implying that the main error in the LDA description is associated with 
local self-interaction effects on B rather than the neglect of the non-analytic discontinuity 
in density functional theory across the band gap. 
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FIG. 1: Electronic band structure of CaBg within the LDA and WDA. 



TABLE I: Theoretical and experimental lattice constant a (A) and internal parameter x for CaLV 





LDA 


WDA 


Experiment 


a 


4.078 


4.131 


4.146 


X 


0.201 


0.206 


0.207 
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